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Abstract. I give a pedagogical review of the derivation for the effective 
lagrangian for nonabelian Debye screening, or hard thermal loops. Following 
Kelly, Liu, Lucchesi, and Manuel, I give the simplest derivation possible, 
using classical kinetic theory. 

The result is valid not just for a thermal, but for an arbitrary initial 
distributions. I use this to study the evolution, at short times, of a gluonic 
"tsunami wave". 

I also suggest how classical kinetic theory may arise at one loop or- 
der. Using the wordline representation of the one loop effective action, I 
follow D'Hoker and Gagne to replace the Wilson line by an integral over 
worldline fermions. A bilinear of these worldline fermions naturally defines 
a nonabelian charge, whose equation of motion is Wong's equation. 



1. Introduction 

In this talk I review our undertanding of plasmas at very high temperatures 
[l]-[9]. While most of the detailed discussion is phrased in terms of a plasma 
near equilibrium at a temperature T, I show that the essential result holds 
not just for a thermal, but for arbitrary, initial distributions of particles. 
The crucial assumption is that the initial density of particles is large. For a 
thermal distribution, this means that the temperature is very high, much 
larger than any mass scale, T ^> m. 

There are two major applications of this work. The first is to the behav- 
ior of QCD in a "deconfined" phase. This is relevant to the central region 
of heavy ion collisions at high energies, as has been explored at the SPS, 
and will be studied at RHIC and LHC. The second is to the electroweak 
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theory, in a phase in which the Higgs phase evaporates. This has become 
of especial interest with regards to baryeogenesis at the electroweak scale. 

All calculations are done using perturbation theory. At zero tempera- 
ture, the normal expansion parameter is the fine structure constant, a = 
g 2 /A7r; in QCD, if one is lucky, a is of order one. For systems at nonzero 
temperature, however, the relevant expansion parameter is not a but the 
coupling constant itself, g. Thus for QCD, if a ~ 1., g ~ 3., and clearly a 
perturbative approach is invalid. It is probably a reasonable approximation 
in the electroweak theory, where g ~ .3. 

My philosophy is that in order to understand any theory, one must 
have some controlled and calculable limit. It certainly should be able to 
give us some sort of picture of the important quasiparticle modes and their 
interactions. In a strong coupling regime the true quasiparticles may behave 
differently, but one must at least know what behavior that they must match 
onto in weak coupling. 

Although it was not clear at first, the basic physical phenomenon which 
I will be studying is that of Debye screening. While it sounds prosaic, even 
to lowest nontrivial order in the coupling constant, for nonabelian gauge 
theories the physics is remarkably involved. For static charges, the screening 
is unremarkable, like that of nonrelativistic systems. It is only when one 
considers the screening of moving charges - that is, the dynamics of the 
system - that a wonderfully rich structure emerges. 

This is suggestive. Numerical simulations of lattice gauge theory have 
clearly demonstrated their power to compute the behavior of gauge theo- 
ries in equilibrium. These numerical simulations have shown that the phase 
structure of QCD depends upon the number and masses of the quark fla- 
vors in a unexpectedly intricate manner. At present, however, lattice gauge 
theory is only efficient for computing the behavior the nature of quantities 
in equilibrium. The present discussion shows that even perturbation theory 
shows that there will be many new suprises in considering truly dynamical 
phenomenon. 

To appreciate the Debye lagrangian one really has to apply it in de- 
tailed calculations. For reasons of sloth, in these proceedings I concentrate 
exclusively on how to derive the effective Debye lagrangian, and not to 
what uses it can be put. The nonabelian Debye lagrangian was first derived 
by brute force, after a laborious and complicated analysis of diagrams in 
perturbation theory [1]. Blaizot and Iancu then derived it by a mean field 
approximation to the Schwinger-Dyson equations, which gives a kind of 
semiclassical kinetic theory [5]. My derivation herein uses classical kinetic 
theory, and copies that of Kelly, Liu, Lucchesi, and Manuel [6]. At one 
point, I use a trick of Brandt, Frenkel, and Taylor [7]. 

I first derive the effective lagrangian in QED, where the gauge principle 



NONABELIAN DEBYE SCREENING 



3 



is elementary, and then that for QCD. The beauty of classical kinetic theory 
is that, after invoking the nonabelian gauge symmetry, the derivation is as 
simple as for the abelian theory. My treatment is meant to be elementary, 
appropriate for those who have never had any exposure to a nonabelian 
gauge theory. This is only possible because of the simplicity of classical 
kinetic theory. 

There is another advantage of classical kinetic theory: it appears to be 
valid not just for thermal, but for arbitrary initial distributions. In sec. 5 I 
use this to study what I term a gluonic "tsunami wave" . 

The apparent peculiarity of classical kinetic theory is that it is necessary 
to introduce a nonabelian charge, Q, and its attendant equation of motion, 
known as Wong's equation [10]. Certainly to most high energy theorists, 
the introduction of Q appears ad hoc and a little peculiar. In the final 
section I outline a possible way in which Q arises directly from the worldline 
formalism. 

2. Classical kinetic theory in QED 

The lagrangian of massless QED is, of course, 

C = H x D a ^ + \F%. (1) 

The covariant derivative and field strength are 

D a = d a - ieA a , F a/3 = d a Ap - d p A a . (2) 
This is invariant under local gauge transformations, 

D a -» ntD a n , v -» nty . (3) 

The transformation of the gauge field is more familiarly written as 

n = , A a ^A a + d a oo . (4) 

For the nonabelian theory, however, it will turn out that thinking about 
covariant derivatives, as in (3), is much more useful than in thinking about 
gauge potentials, as in (4). This difference isn't apparent in the abelian 
theory. 

We then need the classical equations of motion. To be logically complete 
I should start from the lagrangian and derive the appropriate equations in 
the limit of high density. A way in which to do this is sketched in sec. 6. For 
now I will take the equations of motion for granted, since they are utterly 
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standard: introducing the position, x a , and the momentum, p a , there is 
just the definition of the momentum, 

and the Lorentz force equation, 

^=eF aP pf>. (6) 

The field is assumed to be massless, so "r" is not really a proper time, just 
an affine parameter which labels the worldline of the particle. 

As typical of classical kinetic theory, I introduce the density of single 
particles, f(x,p). This is assumed to satisfy the Boltzmann equation in the 
collisionless approximation, 

^f(x, P ) = 0. (7) 

Viewing f(x,p) as a density in phase space, this is a type of Liouville 
equation. 

The only other quantity required is the current. In the classical ap- 
proximation this is merely the product of the charge, the momentum, and 
the density in phase space, integrated over the distribution in momentum 
space: 

f{x) = J d 4 pep a f(x,p). (8) 

Beginning with this expression, one sees how integrals over the distribution 
in momentum space naturally arise. Also, notice that in QED, everything 
— x, p, f(x,p), and j a — are all gauge invariant. 

These elementary equations, (5) - (8), are all that we need to solve for 
the effective action. Using the chain rule, the Boltzmann equation equals 

df fdx a d .d^_d_\ , , . 



dr V dr dx a dr dp a 
Plugging in the equations of motion, 



dr 1 \dx« ap dpP 

which is known as the Vlasov equation. 

Now expand the complete distribution function about some inital value, 

f°, 

f = f° + f + ... (11) 
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The system must start out electrically neutral, with f° the same for positrons 
as for electrons, so initially there is no current. A current is induced by fluc- 
tuations, f 1 . Solve the Boltzmann equation to first order in f 1 : 

df° 

p . d f = ep a Fa(3 J-. (12) 

Without being too careful about what it means, we rather sloppily intro- 
duce the nonlocal operator 1/p ■ d to solve (12) for f 1 , and obtain the 
induced current: 



r-jt^J-rf^w (13) 



Now I follow Brandt, Frenkel, and Taylor [7], and integrate d/dp 7 by parts. 
Remember that the momentum p is just a parameter of the initial (classical) 
distribution, so the field strength tensor is completely independent of it. It 
is then easy to integrate by parts, 

There are three terms from this derivative, which equal 

j =~ e J dp [t^-- 9 WW) Ffhf • ( 5) 

Of the three terms, that ~ 5^ drops out, because it is contracted with the 
antisymmetric field strength, F^. 

Now, implicitly any current defines a lagrangian density through the 
relation 

r = . (i6) 

Thus we can need to find the "Debye" lagrangian, Coebye, which generates 
the current of (15). This is easy: since Cpebye is gauge invariant, it is natural 
to try to use the field strength tensors, ~ F a p, since they are automatically 
gauge invariant. Without too much effort one can see that the result is 

C-Debye = y Jd'p {f^^^F^ f°(p) . (17) 
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3. Classical kinetic theory in QCD 

I defer any discussion of (17) to show that once one sets up the classical 
kinetic theory in a properly gauge invariant fashion, the nonabelian Debye 
lagrangian follows immediately from the abelian. 

So let me start by reviewing what nonabelian gauge invariance is. The 
lagrangian for a massless quark coupled to a SU (N) gauge field looks just 
like that for QED, 

C = iP 7 a D a ^ + \tr(G 2 a/3 ), (18) 

except that now the gauge potential is an SU (N) matrix, A a = A%t a , where 
I normalize that SU(N) matrices t a as tr(t a t b ) = 5 ab /2, with the indices 
a,b = l...(N 2 - 1). 

Normally, textbooks present nonabelian gauge symmetry as a direct 
generalization of the abelian symmetry, giving the transformation of the 
gauge potential, etc. This is a completely confusing way of viewing non- 
abelian gauge invariance. Instead of considering the gauge potential, it is 
much easier to concentrate on the covariant derivative, since under a local 
gauge transformation, it transforms homogeneously: 

D a = d a - igA a -» &D a n , = 1 . (19) 

The transformation of the gauge potential A a can be worked out from this; 
its transformation is inhomogeneous so inelegant. Given (19), however, you 
don't need to worry about the transformation of A a . For instance, if the 
quark fields transforms as in QED, like tp — > Stfif), then obviously the quark 
part of the lagrangian is gauge invariant. 

The field strenth is constructed from the commutator of two covariant 
derivatives. Since D a transforms homogeneously, so does its commutator: 

G al 3 = ^-[D a , D ] - rtG a(S n . (20) 

This is not obvious from the explicit expression for the field strength tensor, 

G af3 = d a Ap - d p A a - ig[A a , A p ] , (21) 

and the transformation of A a . 

The moral of the story is that nonabelian gauge invariance can be easily 
ensured just by sticking to covariant derivatives everywhere. The derivation 
of the nonabelian Debye lagrangian below provides an striking illustration 
of this. 

I now turn to the classical equations of motion for a nonabelian particle. 
As usual the position and momentum are related as p a = dx a /dr. The 
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generalization of the Lorentz force equation of (6) is not trivial, though, 
since while x and p are gauge invariant, in a nonabelian theory (unlike the 
abelian case) the field strenth tensor G a p is not. To form a gauge invariant 
equation, it is necessary to introduce a matrix valued charge, Q: 

^ = 2gtr(QG a p)pP . (22) 

This is trivially gauge invariant if, like G a /3, Q transforms homogeneously 
under a local gauge transformation: 

Q -► ^Qti . (23) 

We now have a different problem - what is the equation of motion for 
Ql This can basically be guessed from gauge and lorentz invariance. We 
certainly want the equation of motion to be gauge covariant, which is easily 
ensured by taking D a Q. This isn't quite right, however, because there are 
then four, instead of one, equation of motion. To get one equation, we con- 
tract D a Q with the obvious vector floating about, which is the momentum. 
This gives us a result first obtained by S. Wong [10]: 

-^D a Q = U. (24) 

The phase space for a nonabelian particle is now x, p, and Q; thus 
the single particle density is in principle a function of all three, f(x,p,Q). 
The classical current is again the product of the charge, g times Q, the 
momentum, and the phase space density: 

j«( x ) = Jd 4 p JdQgQ P a f(x,p,Q). (25) 

I note that the equation of motion for the gauge field is 

D a G^=f. (26) 

Since G a/3 is a commutator of D's, D a DpG a P = 0, so the current is covari- 
antly conserved, 

D a j a = . (27) 

This is only true if Q satisfies Wong's equation. 

Taking these equations for granted, solving for the nonabelian Debye 
action is no harder than for the abelian. The collisionless Boltzmann equa- 
tion, is unchanged: 

J-/(x,p,Q) = 0. (28) 
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The basic point is that because x, p, and f(x,p,Q) are gauge invariant, 
while Q, j a , and the equations of motion are gauge covariant, then the 
resulting Coebye must be gauge invariant. 

Now applying the chain rule to Boltzmann's equation gives a Vlasov 
equation with three terms: 

Tr f = *■ (SS " 2gtr W + 29tr ( lA - Q] lk)) f = ~ <29) 

The last term is from color precession of Q, and is where all of the compli- 
cations of the nonabelian theory reside. 

Suppose we only wish to compute to lowest order in the gauge potential, 
however; this is linear in A a for the current, or quadratic in A a for the 
lagrangian. Then the nonabelian theory is no more complicated than the 
abelian. Expand the distribution function about some initial value, where 
the initial distribution is assumed to depend only upon momentum: 

f(x,p, Q) = f°(p) + f\x,p,Q) + ... . (30) 

Since /o is colorless, to lowest order in A we can drop the term for color 
precession in Boltzmann's equation. In this case the equations abelianize. 
The only feature of the Q's which is needed are the Casimir's: 

J dQ Q a Q b = C 5 ab . (31) 

For gluons in the adjoint representation, C = N, while for quarks in the 
fundamental representation, C = \. 

Thus to ~ A 2 , the nonabelian debye action is a sum over abelian actions: 

4V = 9 2 jd'p tr (^zf^^) E CrfAv) ■ (32) 

Here the field strength is only the abelian part, linear in A a : F a p = d a Ap — 
dpA a . 

(2) 

It is then easy to make C Dehye gauge invariant: replace the abelian part 
of the field strength by the complete field strength, and replace the ordinary 
by the covariant derivative, to obtain: 

£>Debye = 9* Jd 4 p tr ^G af3 ^ G^j ^ C rfr(p) ■ ( 33 ) 

The sum is over the representations of all charged fields. Since everything 
transforms homogeneously under a local gauge transformation, it is obvi- 
ously gauge invariant. The crucial question is whether the generalization is 
unique. On this point I have to defer to the literature [l]-[9]. 
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In hot QED, (17) was first derived almost 40 years ago by Silin [2], 
using classical kinetic theory. For an abelian theory, this is the complete 
result; jC-Debye is quadratic in the gauge potential, and only contributes 
to the photon self energy, which is how Silin wrote it. In constrast, in a 
nonabelian theory, because of the A a which lurks in the covariant derivative 
in l/(p ■ D), there are terms of arbitrary order in the gauge potential. This 
means that C Debye contributes not just to the gluon self energy, but to 
couplings between three, four, or any number of gluons. As a practical 
matter, this complicates using the nonabelian Debye lagrangian in real 
complications, although methods have been developed to deal with this [1]. 

The nonabelian Debye lagrangian was first derived by Taylor and Wong 
[3]. There are several, equivalent ways of writing it; that in (17) is handy 
because it is manifestly gauge invariant [4]. 



4. Thermal distribution 

To understand the significance of the nonabelian Debye lagrangian, consider 
a thermal distribution at a temperature T: 

f° = ^) >< ^T~i > ( 34 ) 

where the — is for the Bose-Einstein distribution of gluons, and the + for 
the Fermi-Dirac distribution of quarks. In this case it is easy to compute 
the nonabelian Debye lagrangian. The integral over the magnitude of p 
generates a factor ~ T 2 , and leaves an angular integral: 

CDebye = \ VT? g J ^ tr ^ ^^ ^ ■ ( 35 ) 

I redefine p — > (i,p), with p a unit spatial vector, p 2 = 1. m g is the gluon 
Debye "mass", 

2 J 9 

There is also a Debye lagrangian for quarks: 

i J <37 » 

where m q = is the quark Debye "mass" , 

2 A^-l^T 2 

m i = ^r — - (38) 



m ' = r + f V' (36) 
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The quark Debye lagrangian is chirally symmetric, like the original gauge 
interaction. 

We can now understand the applicability of the Debye lagrangians. For 
a massless gas at nonzero temperature, the typical momentum is of or- 
der the temperature, which is termed "hard" [1]. The nonabelian Debye 
lagrangian introduces "soft" momenta, on the order of the Debye masses, 
m g ~ niq ~ gT. To understand the relative importance of fields, we use 
power counting, taking the gauge field A a ~ T. Thus for hard momenta, 
the original lagrangian C ~ T 4 , while the nonabelian Debye term is just 
part of the perturbative corrections, Cmbye ~ Q 2 T A . In contrast, for soft 
momenta, both the original and the nonabelian Debye lagrangians are of 
the same order, C ~ Cuebye ~ Q 2 T A . Thus for soft momenta, since Coebye 
is as big as C, it is necessary to include both in an effective lagrangian. 

Diagramatically, one can show that the nonabelian Debye lagrangian 
only receives contributions from one loop diagrams in which the loop mo- 
mentum are hard. This suggests the term "hard thermal loops", which is 
commonly used in the literature. In this talk I adopt instead the more 
generic phrase of nonabelian Debye lagrangian, which sounds less techni- 
cal. I emphasize that while the physics subsumes Debye screening, there is 
much more going on than just that. 

Indeed, consider the textbook example of hot QED. While the photon 
self energy is treated in the textbooks [2], that of hot fermions is not. The 
Debye lagrangian for fermions in hot QED is very similar to that of quarks 
in hot QCD, (37); except for a change in the Debye quark mass, m q , the 
only other change is to use the abelian covariant derivative. As noted, in 
hot QED the abelian Debye lagrangian just contributes to the photon self 
energy, nothing more. But the fermion Debye lagrangian is as complicated 
in hot QED as in hot QCD: the nonlocal factor of l/p-D generates not just 
a contribution to the fermion self energy, but as well an infinite series of 
couplings between a fermion anti-fermion pair and any number of photons. 
Thus the textbook treatment of hot QED is seriously incomplete. There is 
a good reason for this: the fermion Debye lagrangian does not appear to be 
easily derivable by the standard form of classical kinetic theory. It has been 
derived in perturbation theory, and from the semiclassical kinetic theory 
of Blaizot and Iancu [1]; perhaps it could be derived using the approach 
outlined in sec. 6. 

One novel aspect of the Debye lagrangians is that they are nonlocal, 
because of the the factors of 1/p ■ d. In a thermal distribution, this gives 
rise to discontinuities for spacelike momenta, which can be understood as 
a relativistic generalization of Landau damping. 

There is a more general lesson. At zero temperature, a standard as- 
sumption in constructing any effective lagrangian is that all terms must 
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be local. The nonabelian Debye lagrangian shows that this is no longer 
true at nonzero temperature, although the nonlocality which enters is of an 
extremely specific form. 

The nonlocality can be made to (apparently) disappear by choosing the 
gauge p • A = 0. In that case, the entire lagrangian collapses to a mass term 
for the gauge field: 

CDebye = 9* f^P tr {A 2 a ) £ C r /°(p) . (39) 
J r 

This form was first noted by Frenkel and Taylor [1], [7]. Elmfors and Hans- 
son have shown how it can be used to provide a direct functional derivation 
of CDebye from the background field method [8]. Similarly, in this gauge 
the quark Debye lagrangian also collapses to just a self energy. For a gen- 
eral distribution this gauge choice can only be imposed after generalizing 
the gauge potential to be a function not just of spacetime, but also of the 
momentum of the initial distribution, p. 

5. Tsunami wave distribution 

In the above I have made no reference to the initial distribution, /°. It may 
be thermal, but need not be; all that is necessary is that the change from 
the initial distribution is small, f 1 <C /°. For a thermal distribution, a small 
perturbation will automatically return to thermal, but this is certainly not 
true for an arbitrary initial distribution. 

For a general distribution, the nonabelian Debye lagrangian is presum- 
ably valid not just at weak coupling, but only for small times; at long times, 
what /° evolves into depends on the detailed dynamics. I now apply these 
results to another distribution, admittedly cooked up. 

In cascade models of heavy ion collisions, the two colliding nuclei are 
modeled by two pancakes, where in each pancake, all partons move in lock- 
step with the same momentum. Now take one pancake away, and enlarge 
the other pancake until it fills up space uniformly. 

This leads to the "tsunami wave" problem: at time t = 0, assume that 
one has a large and spatially constant density of particles, all moving to- 
gether with the same momentum. Into what state does this evolve at infinite 
time? 

The particles must be bosons, since I have assumed that the initial 
density is high. If the particles have nonzero mass, then one can transform 
into their rest frame, in which case the particles are just a condensate (of 
some sort) at zero momentum. So assume that they are massless, and move 
on the light cone. 

The initial state is a system with nonzero energy and momentum den- 
sity. Thus one natural guess for the state at infinite time is that of a boosted 
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thermal distribution, since that also has nonzero energy and momentum 
density. That the system thermalizes is not obvious; assuming that it does, 
the relevant question then is, over what time scales? 

This is a very difficult question which I could not attempt to solve 
analytically. For small times, however, a perturbative analysis should be a 
reasonable approximation. I define the tsunami wave distribution as: 

f°(p)~ Po 5(p 2 ) \ P 0\SHP-Pb) ■ (40) 

Here po is a parameter proportional to the density, and po = (*|Po|>R>)> 
Pq = 0, is the momentum of the particles in the tsunami wave. 

We now have a problem in which there is a preferred four vector, pq. 
For this problem the obvious choice of gauge is po ■ A = 0. Consider the 
gluon self energy: it is a function of the four momentum, k, and also of this 
vector pq. In general four functions enter into the self energy: 

n^(fc) = u t 5 a ? + u e k a k? + n 3 (k a pP + P %kP) + n 4 P ^ . (4i) 

This is analogous to the situation at nonzero temperature, where the rest 
frame of the thermal bath provides a preferred four vector. For the thermal 
case, the preferred vector is timelike; here it is null. Also as at nonzero 
temperature, the four functions in II Q/3 are related by a Ward identity, but 
that doesn't matter here. 

We can then read off the result from (39). Ignoring inessential constants, 

n t ~+g 2 Np , n £ = n 3 = n 4 = o . (42) 

This is not the whole story, however; it is still necessary to work out the 
effective propagator, including this self energy. This is a straightforward if 
tedious exericse. In po ■ A = gauge, the result is 

A af3 1 ( xaf3 Po^ + k a p^ k 2 PqP q \ 

A = A^TTb \ p^k W^n e {Ut + Ue) WW) ■ 

(43) 

Thus we see that in a tsunami wave, the two, spatially transverse modes 
which one expects for a spin-one field are screened. If IL? were nonzero, it 
would represent a collective mode, analogous to the plasmon of a thermal 
distribution. Because Ilg = 0, there is no plasmon for this type of tsunami 
wave. 

Thus at short times, interactions between gluons in a dense tsunami 
act to screen the usual transverse modes of the gluon. What happens at 
longer times is an open question. This can be studied numerically in a scalar 
theory in the limit of a large number of components [11]. 
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6. Where does Q come from? 

While classical kinetic theory is often used in nonrelativistic systems, its use 
in a nonabelian gauge theory is unfamiliar. What are the original degrees 
of freedom in the gauge theory which conspire, in the classical limit, to 
generate the nonabelian charge Ql 

Let me start with something familiar, the lagrangian for a charged scalar 
field, 4>, in the presence of a background gauge field: 

C = 4> ] {-D 2 )4> . (44) 

Integrating over <\> gives the effective action 

S eff = trlog (-D 2 ) . (45) 



Now I use the usual trick of Feynman and Schwinger, to turn this S e ff 
into the path integral for a particle. First introduce a parameter, r, which 
will become like a proper time: 



dr 

>eff 



Going to momentum space, 

The remaining trace is only for the SU(N) matrix A a . This is then con- 
verting into a sum over paths: 

Seff ~ J^^J VxVp tr (e~ dr ' £(p) ) , (48) 

C(p) = -ip-± + ^(p-gA) 2 . (49) 
where x a = dx a /dr 1 . The integral over the momentum is trivial, 

Seff ~ J^Jvxtr (e~ fo c -^j , (50) 

x 2 

£ A = — -igA- x , (51) 

subject to the boundary conditions x a (r) = +x a (0) 

This form of the one loop effective action has been studied by many 
authors [12]. This representation has been shown to be a powerful way 
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of reorganizing perturbation theory. From (50), this is just a matter of 
expanding the exponent in powers of the gauge field. 

In a more general context, however, (50) is incomplete. In (51), the 
"lagrangian" C is not really that at all: it is a matrix, as in (50) there is 
still a color trace left to do. What we want instead is a form where the 
lagrangian is a true scalar in color space. 

This can be done by replacing the Wilson line by an integral over world- 
line fermions, A(r) [13, 14]: 

tr V e i9 Io A - dx =(-Yj2 f V\V\^ e ^(AtA+iv/2-i) g -/>'£ A (52) 

C x = At (± ■ D) A . (53) 

The A(t) lie in the fundamental representation of SU (N) color, and satisfy 
antiperiodic boundary conditions, ip(r) = —ip(0). Path ordering is denoted 
by V, and is automatic in the path integral. 

It has been known for some time that worldline fermions are required 
to describe finite dimensional representations of nonabelian charge [13]. To 
understand (52), notice that in the absence of a gauge field, the propagator 
for a fermion in one dimension is a step function: the solution to d T A\ = 
S(t) is Aa ~ 9(t — t'). Then, since propagation is in one dimension, the 
complete propagator is naturally an exponential, where the step function 
provides the path ordering. 

The complete form in (52), including the sum over k, was derived by 
D'Hoker and Gagne, eq. (5.2) of [14]. The sum over k is required in order 
to project upon states with occupation number one. 

The integral over worldline fermions then replaces the color sum, and 
gives us a true lagrangian. We can then use this to go back, reintroduce 
the momentum conjugate to the position x (this is not necessary for the 
worldline fermions, since ip^ is already conjugate to ijj), to obtain 

Seff ~ r*Lj VxVpVWtf J2 e ik ^ x+N ' 2 ~V flY e'Io d -' c , 

^~ k=l ^~ 

(54) 

C = -ip ■ x - A^A + ^(p-2g tr(QA)) 2 . (55) 
The nonabelian charge Q is just 

Q a = \U a X . (56) 



The lagrangian C of (55) is similar to the original C(p) of (49), but now the 
color trace downstairs is replaced by an integral over the worldline fermions. 
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Moreover, the nonabelian charge now enters in a most natural manner, as 
a bilinear in the worldline fermions. 

It is easy to check that C gives the correct equations of motion for a 
classical, nonabelian particle. For example, evidently the equation of motion 
for the worldline fermion, x ■ DX = 0, gives Wong's equation for Q, (24). 

A lagrangian very similar to C was proposed by Brandt, Frenkel, and 
Taylor [7] to generate the nonabelian Debye lagrangian. (They used word- 
line scalars instead of fermions, but this difference doesn't matter in the 
classical limit [13].) My contribution is to note that the introduction of the 
worldline fermions is not just a trick to get the correct equations of motion, 
but is a systematic approximation to the correct effective action. 

Given this lagrangian, we then adopt a classical approximation. Instead 
of an integral over the nonabelian charge Q, one should integrate over 
the worldline fermions. Since all we needed before was the Casimir of the 
representation, this doesn't matter in the classical limit. 

Other approximation schemes have been developed to analyze world- 
line path integrals at nonzero temperature. For instance, one can sum over 
paths which wind around in the imaginary time direction [15]. The result- 
ing expressions are not especially simple, though; classical kinetic theory 
appears to be more useful. 

However, this begs the question of how the collisionless Boltzmann equa- 
tion arises. I conclude with a suggestion. In the absence of interactions, at 
nonzero temperature the propagator is directly proportional to the statis- 
tical distribution function: 

tr (An) ~ n W • ( 57 ) 



-a 2 . 

In the presence of a background gauge field, the propagator can be written 
as a sum over paths: 

tr(-^j ~ j~ dr J VxVpVXDtf £ e ik ^ x+N ' 2 ~V e -J>' £ . 

(58) 

From this, define the "single particle" density, /, as 

tT (^2 ) ~ jo ^ l VXT>PVXDX] H X > P ' ^ ' ( 59 ) 

Perhaps in the classical approximation, the dominant term is that where 
/ is stationary with respect to r; this would then give the collisionless 
Boltzmann equation. The expression for the current, (25), follows by dif- 
ferentiation with respect to the gauge field. At two loop order, the effective 
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action is not given by (45), and is how collisions enter into the Boltzmann 
equation. 
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